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I. INTRODUCTION 

Recent research efforts [1-3] have confirmed the ability of the "parabolized" 
Navier-Stokes (PNS) codes to predict accurately and rapidly the aerothermodynamics 
of the actual Space Shuttle Orbiter up to an axial station that is 50% of the overall 
length. This corresponds to a location where the strake ends and the 45° swept wing 
begins. In the past, the geometry of the Orbiter usually has been modified [4-7] 
such that a solution over the complete body was possible. These modifications take 
the form of changing the sweep angle of the wing, removing the canopy, and altering 
the lee of the body so that the cross section is composed of two elliptical shapes. 

In one instance [8], an inviscid solution was obtained for the complete Orbiter. In 
order to perform this calculation, a "fix" had to be used in the vicinity of the bow- 
shock/wing-shock interaction region. Physically, what occurs is a region of embedded 
subsonic flow owing to the bow-shock/wing-shock interaction which causes the result- 
ing coalescing shock wave to be more normal to the flow for a short strearawise dis- 
tance. Since the above inviscid procedure was a marching code, it failed when the 
Mach number m the marching direction became subsonic. 

A viscous numerical procedure is described, to compute the flow over the Shuttle. 
Results are presented that demonstrate the capability of the method. Obtainment of 
these results requires the use of two computer codes. A PNS code [9,10] is used to 
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obtain the solution up to the bow-shock/wing-shock interaction region, and an unsteady 
continuation code is used for the region after the shock interaction. The unsteady 
Navier-Stokes code [11] is also used to obtain the blunt -body starting solution. Only 
results from the marching code will be presented. For the flow conditions calculated, 
that is. Mo, = 7.9, a = 25°, T wa ^^ = 540°R, Re^ = 60,728/in., laminar or turbulent, 
the PNS code has been marched up to an X/L =0.7 which is where the bow-shock/wing- 
shock interaction region occurs. In this work, L refers to the length of the 
vehicle. 

II. COMPUTATIONAL TECHNIQUE 

The PNS equations are obtained from the complete Navier-Stokes equations by 
neglecting the unsteady terms and the streamwise viscous derivative terms. The com- 
plete details of all the terms and derivations can be found in Ref. 12. 

In the present formulation, E, (the marching direction in computational space) is 
a function of x only (axis-normal marching planes) . The governing equations are 
hyperbolic-parabolic in this ^-direction if the inviscid part of the flow field is 
supersonic, if there is no streamwise (axial) separation, and if the pressure gra- 
dient in the viscous region near the wall is treated correctly. However, the system 
of equations still allows for the separation in the crossflow plane (n-4). 



The present PNS code uses the Beam-Warming implicit algorithm to update the 
interior of the region and characteristic, implicit, spatially second-order-accurate 
boundary conditions at the outermost shock wave. An elliptic grid generator of the 
type developed by Steger and Sorenson [13] and further specialized to wing bodies by 
Rai et al. [10] is used to generate the grid for the calculations. 

If the conditions in a particular region are such that the marching procedure 
is invalidated, the unsteady Navier-Stokes (UNS) code is used for these regions. In 
calculating the flow over the Space Shuttle, one such region occurs m the vicinity 
of the bow-shock/wing-shock interaction (a pocket of subsonic flow is encountered) . 

The complete details of the UNS code can be found in Ref. 11. The UNS code is 
extremely versatile and relatively easy to use. It uses either a Beam-Warming 
implicit algorithm or a hybrid scheme due to Rizk and Chaussee [11]. The outer shock 
wave is either fitted or captured. Usually the initial guess is furnished by the PNS 
code, which is modified in some manner to march through regions where it would not 
march before. This procedure is acceptable, since the unsteady code takes this rea- 
sonable guess and iterates in time until a steady-state solution is obtained. 

The domain of this unsteady calculation encompasses the subsonic flow. The 
outflow boundary consists entirely of supersonic axial flow m the mviscid part of 
the flow field. This permits the PNS code, which is more efficient, to continue 
marching from this point. 

III. RESULTS 

Numerical results have been obtained for the following wind-tunnel conditions: 

M^, = 7.9, a = 25°, T wa ^^ = 540°R, Re^ = 60,728/in. turbulent flow. For this calcula- 
tion, the Shuttle surface coordinates were obtained from Rockwell-International 
Corporation. The current geometry consists of the complete Shuttle; the canopy, 

OMS pods, and the vertical stabilizer are included. 

The three-dimensional blunt-body code originally developed by Kutler et al. [14] 
was used to obtain the blunt -nose solution which creates the necessary starting 
planes for the PNS code at X/L = 0.0522. This solution was then marched downstream 
using the elliptic grid generator to construct the grid between the body and the 
fitted outermost shock wave. The grid consisted of either 61 or 121 points in the 
meridional direction and 45 geometrically stretched radial points. An example of the 
grid at an X/L =0.66 is shown in Fig. 1. The outermost grid line is the bow wave, 
which is fitted using an implicit technique. 

The pressure contours m the region of the canopy are presented m Figs. 2 and 3. 
In Fig. 2, the contours on the lee pitch plane of symmetry between X/L = 0.067 and 
X/L =0.4 are shown. In the canopy region, the coalescence of contours details the 
canopy shock wave followed by an expansion wave on the lee of the Shuttle. The con- 
tours for a cross section at an X/L =0.2 are presented m Fig. 3. The canopy 
shock is once again viewed at the point on the lee where the pressure contours 



coalesce. The expansions which are visible on the windward are due to discontinui- 
ties in geometry. 

The Mach number contours at an X/L = 0.66 are presented in Fig. 4. The main 
features are the wing shock and the crossflow shocks on the wing and upper body, 
respectively. These are denoted by the coalescence of the Mach contours. 

In Fig. 5, the crossflow velocity vectors are presented at an X/L = 0.66. Two 
interesting features seen in this figure are the recirculation region in the wing- 
body juncture and the lee vortex. 

The density contours in the vicinity of the wing tip at an X/L = 0.667 are 
shown m Fig. 6. The wing shock and the bow shock have interacted as characterized 
by the bulge in the outer boundary. This is due to the wing shock becoming the outer- 
most surface, with the bow shock being captured. The bow shock appears as the coales- 
cence of the density contours near the outer surface. 

By numerically simulating oil flow on the surface of the vehicle, as m Fig. 7, 
many interesting features are observed. The lines of separation on both the strake- 
wing and the lee of the body are evident by the coalescence of the numerical oil flow. 
The reattachment line is visible on the Shuttle as a series of oil-flow lines diverg- 
ing toward the separation lines . 

The computer-generated particle paths of Fig. 8 exhibit the same trends in the 
flow field that are visible on the Shuttle surface via the oil flow. Specific fea- 
tures are the vortices on the lee which are due to the strake-wing. At this angle of 
attack, the vortices that are generated on the wing impact on the OMS pod. 

IV. SUMMARY 

A procedure has been presented for calculating the flow over vehicles that have 
embedded regions of subsonic flow in the inviscid part of the flow field. A PNS 
marching code is used to obtain the solution up to the bow-shock/wing-shock interac- 
tion region. In this interaction region, the UNS code can be employed since the 
region has a pocket of subsonic flow. Currently, only the results for the marching 
code up to an X/L = 0.667 are included. In the future, the results for the bow- 
shock/wmg-shock region will be available. 
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Fig. 1 Elliptic grid at X/.L *» 0.66. 


Fig. 2 Pressure contours in the lee pitch 
plane of symmetry. 










1 Report No 2 Government Accession No 3 Recipient's Catalog No 

NASA TM 85977 


4 Title and Subtitle 


VISCOUS COMPUTATION OF A SPACE SHUTTLE FLOW FIELD 

6 Performing Organization Code 

A-9799 

Denny S. Chaussee, Yehia M. Rizk* and 

8 Performing Organization Report No 

Pieter G. Bumng 

10 Work Unit No 

9 Performing Organization Name and Address 

T-6458 

Ames Research Center 
Moffett Field, CA 94035 

11 Contract or Grant No 


13 Type of Report and Period Covered 

12 Sponsoring Agency Name and Address 

Technical Memorandum 

National Aeronautics and Space Administration 
Washington DC, 20546 

14 Sponsoring Agency Code 

505-31-01 

15 Supplementary Notes 

^Informatics General Corporation, Palo Alto, CA 

Point of contact. Denny S. Chaussee, Ames Research Center, MS 292A-14 
Moffett Field, CA 94035, (415) 965-6742 or FTS 448-6742 


16 Abstract 


A procedure is presented, as well as some results, to calculate the flow 
over the wxnged Orbiter. This necessitates the use of two computer codes. 

A parabolized marching Navier-Stokes code is used to obtain the solution up 
to the bow shock-wing shock interaction region and for the region after 
the interaction. An unsteady Navier-Stokes code is to be used in the region 
of the shock interaction. Only results for the marching code are presented. 
For the flow conditions calculated, M^ = 7.9, a = 25°, T W3 ;q = 540°R, 

Re L = 60728 per inch, laminar or turbulent, the PNS code has been marched up 
to an X/L = 0.7 which is where the bow shock-wing shock interaction region 
occurs . 


17 Key Words (Suggested bv Author Isl ) 18 Distribution Statement 

Computational fluid dynamics 

Supersonic flow Unlimited 

Shuttle 

Three-dimensional viscous flow Subject category: 34 


19 Security Classif (of this report) 

20 Security Classif (of this page) 

21 No of Pages 

22 Price* 

Unclassified 

. 

Unclassified 

8 

AO 2 


*For sale by the National Technical Information Service Springfield Virginia 22161 





















End of Document 



